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Abstract 

This paper presents two main results. The first result indicates that in materials with broadly 
distributed microscopic heterogeneities, the fracture strength distribution corresponding to 
the peak load of the material response does not follow the commonly used Weibull and 
(modified) Gumbel distributions. Instead, a lognonnal distribution describes more ade- 
quately the fracture strengths corresponding to the peak load of the response. Lognormal 
distribution arises naturally as a consequence of multiplicative nature of large number of 
random distributions representing the stress scale factors necessary to break the subsequent 
"primary" bond (by definition, an increase in applied stress is required to break a "pri- 
mary" bond) leading up to the peak load. Numerical simulations based on two-dimensional 
triangular and diamond lattice topologies with increasing system sizes substantiate that a 
lognormal distribution represents an excellent fit for the fracture strength distribution at 
the peak load. The second significant result of the present study is that, in materials with 
broadly distributed microscopic heterogeneities, the mean fracture strength of the lattice 

system behaves as /// = r^^p + f > an d scales as p,j rs r^~^yj? as the lattice system 
size, L, approaches infinity. 
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1 Introduction 



It is well known that fracture properties and breakdown behavior are very sensi- 
tive to the microstructural details of the material JH]. In ductile materials, grain 
boundaries and second-phase particles form the important microstructural details 
for fracture and damage evolution. In quasi-brittle materials such as ceramics, the 
fracture properties are usually dominated by the size and spatial distribution of 
microcracks, which are often the artifacts of material processing techniques. The 
fracture strength distribution of a quasi-brittle material is significantly influenced 
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by the distribution of microcracks. Probabilistic life design methodologies, "use- 
ful" service life predictive models, and failure risk analysis of a structural material 
component utilize fracture strength distributions in assessing the safety and relia- 
bility of the material component. 

Traditionally, Weibull and (modified) Gumbel distributions based on "weakest- 
link" approach have been widely used to describe the strength of brittle materials. 
These distributions naturally arise from the extreme- value statistics if one assumes 
the following conditions for the distribution of defect cluster sizes in a randomly 
diluted network [|2J]: (1) defect clusters are independent of each other, i.e., they do 
not interact with one another; (2) system failure is governed by the "weakest-link" 
hypothesis, and (3) there exists a critical defect cluster size below which the system 
does not fail, and it is possible to relate the critical size of a defect cluster to the 
material strength. Moreover, if the defect cluster size distribution is described by a 
power-law, then the fracture strength obeys Weibull distribution, whereas an expo- 
nential defect cluster size distribution leads to the Gumbel distribution for fracture 
strengths. However, in heterogeneous materials with broad distribution of disor- 
der, Weibull and Gumbel distributions may not adequately represent the fracture 
strengths corresponding to the peak load response. There are two main reasons be- 
hind this inadequacy. 

First, in the "weakest-link" hypothesis, the fracture strength of a randomly diluted 
lattice system is determined by the presence of few critical defect clusters, and is de- 
fined as the stress required for breaking the very first "weakest-link" in the system. 
In materials with broad disorder, the breaking of the very first bond ("weakest- 
link") does not usually lead to the entire system failure, and hence the fracture 
strength distribution based on the very first bond failure may not be applicable. In 
addition, the evolution of initial distribution of defect clusters, i.e., subsequent bond 
failure is controlled not only by "weakest" bonds with smallest thresholds but also 
by the stress concentration and shielding effects around the defect clusters. Con- 
sequently, at the peak load, the defect cluster size distribution that evolved under 
the applied stress field may be quite different from the initial defect cluster size 
distribution. Thus, unless the defect cluster size distributions at the peak load and 
at the very first bond failure have the same form, it is highly unlikely that the frac- 
ture strength distributions corresponding to the peak load and the first bond failure 
would follow the same distribution. 

Secondly, if the defect cluster size distribution at the peak load exhibits multi-modal 
behavior then the corresponding fracture strength distribution is not adequately de- 
scribed by unimodal distributions such as Weibull and Gumbel distributions. In the 
presence of experimentally observed multiple defect populations, the tail cannot be 
fitted by a power law with a single power law exponent due to the superposition 
of different defect cluster size distributions. These multi-modal defect cluster size 
populations result in multi-modal strength distributions ||3j.|4|,|5|,|6|,|2|]. Experimen- 
tally, bimodal Weibull strength distributions were observed in various brittle ma- 
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terials including carbon [7] and silicon carbide fibers [6], and for certain ceramics 
|0]. The effect of multiple flaw (defect) populations on fracture strength distribu- 
tion is also examined experimentally in certain ceramic materials il EH EH EH]- In 
general at the peak load, the broken bond cluster distribution follows generalized 
Gamma-distribution and in this case, neither the Weibull nor the Gumbel distribu- 
tions fit the fracture strength distribution accurately 

This paper presents two main results. The first result indicates that, in the case of 
materials with broadly distributed heterogeneities, a lognormal distribution repre- 
sents the fracture strength of the macroscopic system more adequately than the 
conventional (modified) Gumbel and Weibull distributions. The second result indi- 
cates that the mean fracture strength decreases with increasing lattice system size, 
L, and scales as ( Lo g L ^ m the thermodynamic limit. 

The paper is organized as follows. Section 2, presents the theoretical derivation for 
the fracture strength distribution of a lattice system. In Section 3, numerical sim- 
ulations using two-dimensional triangular and diamond (square lattice inclined at 
45 degrees between the bus bars) lattice networks are presented. Section 4 presents 
the validation results and a comparison study between the lognormal distribution 
for fracture strengths derived in this study versus the traditionally used Weibull and 
Gumbel distributions. 



2 Fracture Strength of a Discrete Lattice System 



Progressive damage evolution leading to failure of disordered quasi-brittle mate- 
rials has been studied extensively using various t ype s of discrete lattice models 
lE4lE5l.E6l.E7l ]. Electrical fuse/breakdown models tl 1 8L fiol l2X^ I2TL I22L [23I1 . central- 
force models|24 , ^HHBSB, bond-bending models T3llH 7and beam- 
type models H151 |33J have been used in combination with disorder either in the 
elastic constants, threshold values or in the random dilution of the bonds to model 
damage evolution in brittle materials. The reader is referred to 1 17, 14, THEi! and 
the references therein for a comprehensive review of modeling quasi-static progres- 
sive damage evolution in brittle materials using discrete lattice networks. Although 
discrete lattice models do not describe the specific behavior of any real material, 
they incorporate the essential ingredients of a breakdown process, namely, the ini- 
tial material disorder and the redistribution of forces due to damage evolution. In 
this respect, any realistic progressive damage evolution model that describes the 
behavior of real materials should be capable of reproducing the behavior of these 
idealized discrete lattice models |34]. 



The essential features of discrete lattice models are disorder, elastic response char- 
acteristics, and a breaking rule for each of the bonds in the lattice. The elastic 
and breaking response characteristics of each bond in the lattice correspond to the 
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mesoscopic response of the material. The elastic response of the individual bonds is 
typically described by electrical fuse models, central-force (spring) models, bond- 
bending spring models, and beam-type models. The quenched disorder in the sys- 
tem is introduced either in the elastic constants, in the threshold values, or by the 
random dilution of bonds. The breaking of a bond occurs irreversibly when the ap- 
plied action (stress or displacement) across the bond exceeds the breaking thresh- 
old. Various types of breaking rules have been adopted in the literature depending 
on the elastic response characteristics of individual bonds. When an external action 
(displacement or force) on the lattice is slowly increased, the individual bonds in 
the lattice will break irreversibly one after another until the system falls apart. It 
is supposed that successive fuse failures leading ultimately to the failure of lattice 
system is similar to the breakdown of quasi-brittle materials. 

In this study, we consider a discrete lattice system in which all the bonds are intact 
at the beginning of the analysis and the damage is accumulated progressively by 
breaking one bond at a time until the entire lattice system falls apart. This approach 
is in contrast with the earlier works, wherein either randomly diluted lattice sys- 
tems close to the percolation threshold were considered or the stress required to 
break the very first bond was defined as the fracture strength of the lattice system. 
By using either of these methodologies, it is possible to analysize very large lattice 
systems of size 1000 x 1000. However, the defect cluster distribution obtained by 
randomly diluting the bonds of a lattice systems is quite different from that ob- 
tained by sequentially breaking one bond at a time. This is because, the evolution 
of initial distribution of defect clusters is controlled not only by "weakest" bonds 
with smallest thresholds but also by the stress concentration effects around the de- 
fect clusters. In this paper, we define the fracture strength of a lattice system as the 
stress corresponding to the peak load of the lattice system response. 

Consider a lattice system with a total number of bonds, N el , subjected to a stress 
controlled loading. Let M = {1, 2, 3, - • • , N e i} denote the set of individual bonds 
in the lattice system. After breaking k number of bonds, let S% denote the set of all 
the k number of bonds that were broken. A broken bond is considered a "primary" 
bond if an increase in the applied stress is necessary to break that bond. Similarly, 
let S% denote the set of bonds that will be broken ("avalanche" bonds) without 
any further increase in the externally applied stress, given the set of broken bonds 
S\. at an applied stress a. Figure 1 illustrates this concept using a typical lattice 
system response. Once a "primary" bond is broken, subsequent breaking of the 
"avalanche" bonds continues until we encounter another "primary" bond, i.e., until 
it is required to increase the applied stress level to break a bond. If we do not 
encounter a "primary" bond, it means that we reached the peak load of the lattice 
system. The dimension of the set S% represents the avalanche size, and is given by 
Ofe. Let <Sfc denote the set of unbroken bonds remaining in the lattice system after 
k number of bonds and the elements of the avalanche set S% have already been 
broken. It should be noted that the set <S£ also includes the set of bonds T>% called 
the dangling or dead ends that do not carry any stress even though they are not 
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broken. The sets S\, S k and S k are mutually disjoint and collectively exhaustive 
such that 



s°ns% = 
si n si = 
s b k n si = 



(i) 



where denotes the null set, and © denotes the additive sum of disjoint sets. As- 
sume that for each non-empty set S k , a set of non-negative stress concentration 
constants C k = {\{S\ © S k ), i e S k } are defined. These constants represent the 
stress concentration coefficients for all the elements of the set of unbroken bonds 
S k when all the elements of the sets S^ and S k are removed from the lattice system. 
Denote the stress concentration coefficient Xi(S k © S k ) to be equal to zero for all 
the bonds {i e T> k }. These stress concentration coefficients C k implicitly define 
the redistribution of stress within the lattice system after all the elements in the set 
S k © S k have been removed. 

Let Mk denote the set of non-negative constants such that M. k = {fa = 1/Aj : 
i E (S% - VI)}. Also, let = {x, t = ^of : m G M k , of ~ P^{o) and i E 
(S k — T> k )} denote the set of modified breaking thresholds (similar to annealed 
disorder) for all the bonds in the set (S k —V k ) after the bonds in the set S k ®S k have 
been broken. In the above description, PQ h (cr) denotes the cumulative probability 
distribution of the breaking thresholds at th state, i.e., before any of the bonds 
in the lattice system are broken. However, as the bonds are broken successively, 
the distribution of breaking thresholds P k h (cr) of the set of elements X k h changes 
gradually from its initial distribution Pq H (<j) of the elements in M . 

The conditional probability that a bond from the set S' k breaks given that the bonds 
in the set <S£ © S k have been broken, i.e., for breaking the (k + a k + l) th bond in 
the lattice system given that (k + a k ) have been broken already, is given by 

/(fe+o fc +l) = Pk h ( a {k+a k +l)) ( 2 ) 



where f/ k+ak+ i> denotes the conditional probability of breaking the (k + a k + l) th 
bond, and cr^ 0fc+1 ) = inf(X k h ). Note that a|£+ afc+1 ) denotes the minimum exter- 
nally applied stress that is necessary to break the (k + a k + l) th bond. The condi- 
tional probabilities for breaking any of the bonds within the set S k , given that all 
the bonds in the set S k are broken is equal to one, i.e., all the bonds in the set S k 
break without any further increase in stress after breaking the k th bond. Hence, 

f(k+j) = 1 Vj = {1,2, ■■■,a k } (3) 
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Thus, the probability of breaking the (k + a k + 1) bond in the lattice system is 

F(k+a k +l) = f(k+a k +l) F(k+a k ) 

= f(k+a k +l) f(k+a k ) f(k+a k -l) ' ' ' f(k+l) F( k ) 

= f(k+a k +l) F( k ) (4) 

Let B denote an ordered set of indices, where each index refers to the sequential 
number of the broken bond for which an increase in the applied stress is required 
to break the bond. Then, this set contains the indices such as k and (k + a k + 1) as 
its elements, and is written in a general way as 

B = {b j+1 = (bj + aj + 1) : b = 0; a = 0; j G {0, 1, 2, • • •}} (5) 

In Eq. (5), a,- refers to the avalanche size after breaking the bj indexed bond. It 
should be noted that the set B maps the sequential number of the "primary" broken 
bond to the sequential number of the broken bond. That is, bj = B(j), where j = 
{1, 2, 3, • • •} refers to the sequential number of the "primary" broken bond and bj 
refers to the sequential number of the broken bond. With this notation in hand, Eq. 
(4) can be written recursively as 



F (k+a k +l) = n f(j) (6) 

Similarly, let A denote an ordered set of externally applied stress values that are 
required to break the "primary" bonds, i.e., 

A={af n : jeB} (7) 



Now, consider the set of scale factors, Q, defined as 

G = K 477 : ''J BU>, bj-i = B{ 3 — 1); gi = 1; j = {2, 3, • • •}} (8) 

Since the elements of the set A depend not only on the initial distribution of break- 
ing thresholds but also on the stress concentration factors around the broken bond 
clusters, these stress levels, <7™ n , at which the "primary" bonds break can be con- 
sidered as independently distributed random variables. Consequently, the elements 
of the set Q are also independently distributed random variables. Using Eqs. (7) and 
(8), the stress required to break the (k + a k + l) th bond, a}%™ ak+ u, can be expressed 
as a product of independently distributed random scale factors, i.e., 

vain I TT 1 rain /c\\ 

°(k+a k +l) = I 11 93 J a l ( 9 ) 
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where o™ m is the stress required to break the first bond. 

In lattice systems with broadly distributed breaking thresholds, the cumulative avalanche 
sizes in the regime up to the peak load of the material response are negligible com- 
pared with the total number of bonds broken up to the peak load. Consequently, 
the dimension (cardinality) of the set Q with independently distributed scale factors 
gj, j E B, is approximately 0(n p ), where n p m 0(L 1,8 ) is the number of broken 
bonds up to the peak load. By virtue of the central limit theorem, which states that 
the product of a large number of independent factors, none of which dominates the 
product, will tend to the lognormal distribution regardless of the distributions of the 
individual factors involved in the product, we have, Prob[cr^™ afc+1 ) < a] ~ LN as 
dim(i3) — > oo, where LN denotes lognormal distribution. Hence, as the number 
of "primary" broken bonds at the peak load increases with increasing system sizes, 
the fracture strength distribution for larger lattice systems tends to be lognormal 
distribution. 

REMARK 1: In the case of narrowly distributed breaking thresholds, breaking of a 
bond significantly influences the subsequent bond breaking process, and the lattice 
system reaches its peak load soon after breaking fewer "primary" bonds. In partic- 
ular, in the weakest-link hypothesis, the fracture strength of the system is identified 
with the breaking of the first bond, and hence the fracture strength distribution is 
based on extreme- value theory. Qualitatively, the assumption of broadly distributed 
breaking thresholds is supposed to distinguish the scenario of large number of "pri- 
mary" broken bonds from the case of fewer "primary" broken bonds before the 
system reaches peak load. It is in this sense that we use the notion of broadly dis- 
tributed heterogeneities. 



3 Numerical Simulations 

In this study, we pursue an electrical equivalence to the mechanical problem 1E4II35I 
Bill . We assume equivalence between electrical current, voltage, and conductance 
in the electrical system and the mechanical stress, strain, and Young's modulus in 
the mechanical system, respectively. Similarly, we also assume an equivalence in 
the breakdown process, that is, equivalence between successive burning of fuses 
leading to the loss of electrical network conductivity and the mechanical break- 
ing of bonds leading to lattice system failure. The main advantage of modeling a 
mechanical problem using an electrical analogy is that the number of degrees of 
freedom in the system is significantly reduced, thereby increasing the system size 
that can be simulated using the same computational power. 

Consider a two-dimensional lattice of size L x L. The model adopted in this work 
is similar to the model III (random threshold) presented in Ref. [25], however, as 



7 



mentioned in Section 2, we start with a fully intact lattice system and break one 
bond at a time until the lattice system falls apart. In this work, numerical simu- 
lations are performed on triangular and diamond lattice topologies with periodic 
boundary conditions in the horizontal direction. The elastic response of each bond 
in the lattice is linear up to an assigned threshold value at which brittle failure of 
the bond occurs. The disorder in the system is introduced by assigning random 
maximum threshold current values, i c , (which is equivalent to the breaking stress 
in mechanical problem) to each of the fuses (bonds) in the lattice, based on an 
assumed probability distribution. The electrical conductance (stiffness in the me- 
chanical problem) is assumed to be the same and equal to unity for all the bonds in 
the lattice. This is justified since the conductance (or stiffness) of a heterogeneous 
solid converges rapidly to its scale independent continuum value. The probability 
distribution of failure thresholds is dependent on the particular type of material 
considered. However, since our focus here is not on modeling any specific type of 
material, but on the generic features of damage evolution in disordered systems, 
we choose, a uniform probability distribution, which is constant between and 
1 . A broad thresholds distribution represents large disorder and exhibits diffusive 
damage (uncorrelated burning of fuses) leading to progressive damage localiza- 
tion, whereas a very narrow thresholds distribution exhibits brittle failure in which 
a single crack propagation causes material failure. 

Periodic boundary conditions are imposed in the horizontal direction to simulate an 
infinite system and a constant voltage difference (displacement or strain) is applied 
between the top and the bottom of lattice system. The simulation is initiated with 
a lattice of intact fuses in which disorder is introduced through random breaking 
thresholds. The voltage V across the lattice system is increased until a fuse (bond 
breaking) burns out. The burning of a fuse occurs whenever the electrical current 
(stress) in the fuse (bond) exceeds the breaking threshold current (stress) value 
of the fuse. The current is redistributed instantaneously after a fuse is burnt. The 
voltage is then gradually increased until a second fuse is burnt, and the process is 
repeated. 

The above choice of redistributing the current after breaking a fuse assumes that the 
current relaxation in the lattice system is much faster than the breaking of a fuse. 
Thus, each time a fuse is burnt, it is necessary to calculate the current redistribution 
in each of the fuses in the lattice. This is very time consuming, especially with in- 
creasing lattice system size. The authors have developed a multiple-rank Cholesky 
updating algorithm for modeling relaxation processes in disordered systems W37I1 . 
In comparison with the most sophisticated Fourier accelerated iterative schemes 
used for modeling lattice breakdown [38, 39], this algorithm significantly reduced 
the computational time required for solving large lattice systems. Figure 2 presents 
the snapshots of damage evolution for the case of a uniformly distributed random 
thresholds model problem in a triangular lattice system of size L = 512. Based on 
these snapshots (Fig. 2(a)-(e)), it is clear that the bond breaking occurs more or less 
randomly until very close to the peak load. Since the response of the lattice system 
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based on the above numerical algorithm corresponds to a specific realization of 
random bond breaking thresholds, an ensemble averaging of the numerical results 
is necessary to obtain a realistic representation of the lattice system response. Table 
1 presents the number of configurations, N con f ig , over which statistical averaging 
is performed for different lattice sizes. 



4 Results and Discussion 



4.1 Distribution of Fracture Strengths 



Conventionally, Weibull and Gumbel distributions are used to fit the fracture strength 
data for brittle materials 0.0]. However, as Weibull mentioned in his pioneering 
paper, the Weibull distribution should be considered as an empirical one on an equal 
footing with other type of distributions 1 4(J . In material science applications, log- 



normal, power law, Gamma, Type-I extreme value, and bimodal distributions are 
also often used for describing the fracture strength distribution Jalfl |& 

In randomly diluted disordered elastic networks, a Gumbel distribution better fits 
the fracture strengths distribution far away from the percolation threshold and a 
Weibull distribution provides a better fit close to the percolation threshold. These 
ideas are generally based on the functional form of probability density of the defect 
clusters that are obtained from broken bond cluster statistics Jl7, la, 18, 1^, 2pll . 



In these randomly diluted disorder problems, the defect cluster size distribution is 
exponential far away from the percolation threshold and follows a power law close 
to the percolation threshold. An exponential defect clusters distribution leads to 
Gumbel distribution for fracture strengths and a power-law distribution of defect 
clusters leads to Weibull form for fracture strengths. 



Duxbury et al 11811 studied the distribution of fracture thresholds in the disordered 
media with randomly diluted bonds. Subsequently, this study was applied to a va- 
riety of lattice models with different local behavior including central-force spring 
models |2(J 13], Born models |41], and bond-bending models lUill . In all these 



studies, fracture strength of the lattice system was defined as the stress required for 
breaking the very first bond in the system. These studies along with the analytical 
investigations based on largest crack size W22I1 concluded that the fracture thresh- 
old decreased with the system size L as a power law of log(L) and the fracture 
threshold distribution is best described by a double exponential (modified Gumbel) 
distribution. However, numerical simulations on intact two-dimensional triangular 
and diamond lattice topologies for heterogeneous materials with a broad disorder 
indicate that the entire lattice network does not fall apart as soon as the first bond 
is broken. In fact, the number of broken bonds at the peak load scales as a power 
law. This suggests that the arguments used for deriving the breaking strengths of 
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the first bond may not be extended to this type of problems. 



In this study, fracture strength of a lattice system is defined as the stress correspond- 
ing to the peak load of the lattice system response. A schematic of a typical lattice 
system response for a given distribution of random bond breaking thresholds is 
shown in Fig. 1. Each realization of the random bond breaking thresholds results in 
a specific fracture strength value. The distribution of these fracture strengths sam- 
pled over an ensemble of configurations is the subject of interest in this section. 
In the following, we investigate the validity of Weibull and Gumbel distributions 
to represent the fracture strength data corresponding to the peak load of the lattice 
response. 

The Gumbel distribution for fracture strengths Of is given by 

P G (a f ) = 1 - exp(-cL d exp(~)) (10) 

a f 

and the Weibull distribution is 

P w (a f ) = 1 - exp(-cL d af) (11) 

where k, 5, c and m are constants, and d denotes the lattice dimension. The validity 
of Gumbel and Weibull distributions to the fracture strength data can be tested by 
rewriting Eqs. (10) and (1 1) as 



(12) 

for the Gumbel distribution and 





- Inc (13) 
for the Weibull distribution. In Eqs. (12) and (13), the variable A is defined as 



A = -ln 



ln(l-P(a f )) 



(14) 



where P(<Jf) is Pai^f) m the case of Gumbel distribution, and P w ((jf) for Weibull 
distribution. Figure 3(a) presents Gumbel fit for the fracture strength distribution for 
triangular lattice network using Eq. (12). Similarly, Fig. 3(b) presents the Weibull 
distribution fit (Eq. (13)) for the fracture strengths obtained using simulations on 
triangular lattices. From these figures, it is clear that fracture strength data obtained 
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for two different lattice sizes does not align onto a single straight line as it should, 
if the data were to follow Eq. (12) or (13). The fracture strength results based on 
diamond lattice network exhibit similar trends indicating that in the case of highly 
disordered materials, neither Gumbel nor Weibull distributions may represent the 
fracture strengths distribution accurately. 

The first main result of this study is that for materials with broadly distributed 
heterogeneities, a lognormal distribution represents the fracture strength of the 
macroscopic system more adequately than previously used (modified) Gumbel and 
Weibull distributions. Lognormal distribution has been used in engineering practice 



to represent the fracture strengths of materials W42L 14311 . The lognormal distribution 
can be understood to have evolved as a consequence of multiplicative nature of 
large number of random distributions representing the individual scale factors nec- 
essary to break the subsequent "primary" bonds leading up to the peak load. The 
precise character of the individual distributions that are multiplied to give the final 
distribution is irrelevant as long as the number of "primary" broken bonds up to the 
peak load is large. Figures 4(a) and 4(b) present the cumulative fracture strength 
versus the standard lognormal variable, £, defined as £ = Ln ( (T /)~ T? > f or triangular 
and diamond lattice networks respectively. In the above description, r] and ( refer 
to the mean and the standard deviation of the logarithm of 07. These figures indi- 
cate that the fracture strength distribution collapses onto a single curve for different 
lattice system sizes, which is an improvement compared to (modified) Gumbel and 
Weibull distributions. A better representation to test the lognormal description for 
fracture strengths is to plot the inverse of the cumulative probability, $ _1 (P(erf)), 
against the standard lognormal variable, £. In the above description, $(-) denotes 
the standard normal probability function. Figures 4(c) and 4(d) present the lognor- 
mal fit for the cumulative fracture strength distributions obtained for triangular and 
diamond lattice networks respectively. From these figures, it is clear that the frac- 
ture strength distribution obtained for different lattice system sizes collapses onto a 
single curve, albeit minute deviation from straight line behavior is evident. We have 
also used the normal distribution to collapse the fracture strength data of triangular 
and diamond lattice systems. Although the data collapse is reasonable, it is not as 
good as that of lognormal distribution. 



4.2 Mean Fracture Strength 



The second main result of this study is concerned with the scaling law for the 
mean fracture strength of the lattice system. Table 1 presents the mean fracture 
strength data for various triangular and diamond lattice system sizes. Cumulative 
fracture strength distributions may be used to derive the scaling form of mean frac- 
ture strength. For example, when the Weibull distribution represents the fracture 
strengths, the mean fracture strength, jif, scales as a power law, i.e., 
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d 



Hf~L ™ (15) 

Similarly, in the case of fracture strengths represented by Gumbel distribution, we 
have 



A- 



k 



(In c — ai) + d In L 

1 1 

= A 1 + B 1 lnL = ln(A 2 L Bl ) (16) 

where a\, Ai, B\, and A 2 = exp(Ai) are constants that are related to the param- 
eters k and c of the Gumbel distribution. Equation (16) is same as the relation 
proposed by [2(1 2^] that shows the effect of sample size L on average breaking 



voltage or failure strength. Figures 5(a) and 5(b) present the mean fracture strength 
fj,f versus the lattice size L based on the form of Eqs. (15) and (16), respectively, for 
diamond lattice system. The results for triangular lattice systems exhibit similar be- 
havior. From these figures, it is evident that mean fracture strength does not follow 
a power-law dependence on the lattice size L. The results presented in Fig. 5(b) in- 
dicate that Eq. (16) may be able to represent the mean fracture strengths reasonably 
well even though the Gumbel distribution is inadequate to represent the cumulative 
fracture strength distribution. The value of the exponent S is approximately equal to 
2.45 and is consistent with the values reported in the literature [20, 26]. The coeffi- 
cient B\ (slope in Fig. 5(b)) is large 20.5) and is, once again, in agreement with 
the divergence behavior expected as the lattice system approaches failure |2jJ 2(J. 



This weak (logarithmic) dependence of mean fracture strength on system size L 
is consistent with the analytical and previously reported numerical results on ran- 
domly diluted networks llHfliflifl^EiEill. 



Alternatively, we have plotted the ensemble averaged peak load, F pea k, versus the 
lattice system size, L, as shown in Fig. 6. From the numerical simulation results 
presented in Fig. 6, it is clear that the ensemble averaged peak load F peak for triangle 
and diamond lattice topologies may be expressed as 



peak ' 



C L a + d (17) 



where C and C\ are constants. Thus, the mean fracture strength, jjt,f, defined as 



F 

fif — -sp-, is given by 



u s = CoL a - x + ^ (18) 

The results shown in Fig. 6 indicate that for both the triangular and diamond lattice 
topologies, the exponent a in Eq. (17) is approximately equal to 0.96. This, in 
turn results in a very small exponent value equal to -0.04 in the first term of the 
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Eq. (18). A very small negative value of the exponent (a — 1) is equivalent to 
a logarithmic correction, i.e., for (1 — a) « 1, L a ~ l ~ {log(L))~^ . Thus, an 
alternative expression for the mean fracture strength may be expressed as 



"'=(W + z (19) 

where and c are constants that are related to the constants C and C\ of Eq. 
(18). This shows that the mean fracture strength of the lattice system decreases 
very slowly with increasing lattice system size, and scales as /ij jj^zyj; f° r ver y 
large lattice systems. 



5 Conclusions 



This paper presents a theoretical investigation supplemented by numerical simula- 
tions to describe the fracture strength distribution of a lattice system. The discrete 
lattice system considered is fully intact at the beginning of the analysis and the 
damage is accumulated progressively by breaking one bond at a time until the en- 
tire lattice system falls apart. The fracture strength of a lattice system is defined 
as the stress corresponding to the peak load of the lattice system response. This is 
in contrast with the earlier works, wherein either randomly diluted lattice systems 
close to the percolation threshold were considered or the stress required to break 
the very first bond was defined as the fracture strength of the lattice system. 

Our study presents two main results. First, for materials with broadly distributed 
heterogeneities, a lognormal distribution represents the fracture strength of the 
macroscopic system more adequately than previously used distributions such as 
(modified) Gumbel and Weibull. The lognormal distribution can be understood to 
have evolved as a consequence of multiplicative nature of large number of random 
distributions representing the stress scale factors necessary to break the subsequent 
"primary" bonds leading up to the peak load. The precise character of these indi- 
vidual distributions that are multiplied to give the fracture strength distribution is 
irrelevant as long as the number of "primary" broken bonds up to the peak load 
is large. Hence, as a consequence of the central limit theorem, the system fracture 
strength probability distribution approaches a lognormal distribution. Numerical 
simulations based on two-dimensional triangular and diamond lattice topologies 
substantiate that a lognormal distribution represents an excellent fit for the fracture 
strength distribution. 

Second, the mean fracture strength of the lattice system behaves as fif = ^* f L ^ + f> 
and scales as fif « n^Ttf as me lattice system size, L, approaches infinity. 
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Table 1 
Peak Load 



L 


N C onfig 


Triangular 


Diamond 


Mean 


Std 


Mean 


Std 


4 


50000 


2.087 


0.4328 


1.710 


0.4248 


8 


50000 


3.148 


0.4545 


2.338 


0.4209 


16 


50000 


5.261 


0.518 


3.650 


0.4485 


24 


50000 


7.336 


0.5822 


4.953 


0.4865 


32 


50000 


9.377 


0.6434 


6.244 


0.5219 


64 


50000 


17.372 


0.8595 


11.301 


0.6611 


128 


12000 


32.82 


1.2952 


21.09 


0.9045 


256 


1200 


62.79 


2.0251 


40.06 


1.380 


512 


200 


120.49(121.13*) 


3.5789 







★ estimated value based on the equation given in the inset of Fig. 1 
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• "Primary" broken bonds 
'Avalanche" broken bonds 



Peak Load 



O 




Displacement 



Fig. 1. A typical lattice system response indicating "primary" and "avalanche" broken 
bonds. 



19 




Fig. 2. Snapshots of damage in a typical triangular lattice system of size L = 512. Number 
of broken bonds at the peak load and at failure are 83995 and 89100, respectively, (a)-(i) 
represent the snapshots of damage after breaking n b number of bonds, (a) n b = 25000 (b) 
n b = 50000 (c) n b = 75000 (d) n b = 80000 (e) n b = 83995 (peak load) (f) n b = 86000 
(g) n b = 87000 (h) n b = 88000 (i) n b = 89100 (failure) 
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0' 1 1 1 1 1 

-0.5 0.5 1 1.5 2 

ln(1M f ) 

(b) 

Fig. 3. Probability distribution fits for fracture strengths at the peak load response in a 
triangular lattice for different lattice system sizes L = {4, 8, 16, 24, 32, 64, 128, 256, 512}. 
(a) Gumbel distribution (b) Weibull distribution 
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Standard Variate = (x-li)/g 



Standard Variate = (x-u)/a 



(a) (b) 




-8 -6 -4 -2 2 4 -8 -6 -4 -2 



(c) (d) 

Fig. 4. Lognormal distribution fit for fracture strengths at the peak load response (a) tri- 
angular lattice (b) diamond lattice. Reparametrized lognormal distributions: (c) triangular 
lattice (d) diamond lattice. For triangular systems, L = {4, 8, 16, 24, 32, 64, 128, 256, 512}, 
and for diamond lattices, L = {4, 8, 16, 24, 32, 64, 128, 256} 



22 



-0.6 r 



-0.8 ■ 



-1 



-1.2 ■ 



-1.4 - 



-1.6 



-1.8 - 



1 1.5 



100r 



90 - 



80 ■ 



70 - 



60 - 



2.5 3 3.5 

Ln(L) 



4.5 



(a) 



5.5 



c? 50 ■ 



40 - 



30 ■ 



20 ■ 

10 - 



Slope = 20.48; intercept = -18.32 



1 1.5 2 2.5 3 3.5 4 4.5 5 5.5 6 

Ln(L) 

(b) 

Fig. 5. Mean fracture strength versus diamond lattice system size, (a) Weibull fit based on 
Equation (15) (b) Modified Gumbel fit based on Equation (16) 
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Fig. 6. Average peak load versus the lattice system size, L 
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